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TECHNICAL  ME^DRANDUM  I36I+ 


THE  PLANE  PROBLEM  OF  THE  FLAPPING  WING*^ 
By  Walter  Birnbaiom 

In  connection  with  my  report  on  the  lifting  vortex  sheet^  which 
forms  the  essential  basis  of  the  following  investigations,  I  shall  show 
how  the  methods  developed  there  are  also  suitable  for  dealing  with  the 
air  forces  for  a  wing  with  a  circulation  variable  with  time.   I  shall, 
in  particular,  develop  the  theory  of  a  propulsive  wing  flapping  up  and 
down  periodically  in  the  manner  of  a  bird's  wing.  I  shall  show  how  the 
lift  and  its  moment  result  as  a  function  of  the  flapping  motion,  what 
thrust  is  attainable,  and  how  high  is  the  degree  of  efficiency  of  this 
flapping  propulsion  unit  if  the  air  friction  is  disregarded.  Finally, 
I  shall  treat  an  interesting  case  of  dynamic  instability  for  a  spring- 
suspended  wing;  this  phenomenon  was  confirmed  by  experiments  at  the 
Gottingen  aerodynamic  test  laboratory.  Professor  Prandtl  gave  me  his 
guidance  concerning  the  present  report,  and  I  want  to  express  here  also 
my  sincere  gratitude  for  the  abundant  stimulation  and  energetic  assist- 
ance he  gave  to  me  at  all  times . 

1.  General  statements.-  The  calculations  refer  to  the  two-dimensional 
problem,  that  is,  to  the  wing  of  infinite  length,  or  bounded  by  plane  side 
walls,  or  to  a  wing  with  so  large  an  aspect  ratio  that  the  boundary  effect 
is  negligible.  For  the  rest,  the  same  assumptions  as  in  the  first  report 
are  valid  concerning  smallness  of  the  air  forces,  slight  camber  of  the 
wing,  and  so  forth.  The  wing  is  ass-umed  to  extend  from  x  =  -l/2  to 
X  =  +3/2  so  that  the  point  x  =  0  becomes  the  center  of  pressiore  of 
a  plane  wing  with  fixed  angle  of  attack.  Let  v,  in  the  direction  of  the 
positive  X-axis,  be  the  air  velocity  at  a  large  distance.  Let  the  posi- 
tive Y-axis  point  downward.  By  7  =  7(x,t),  I  denote  the  density  of  the 
lifting  vortices  (now  a  function  of  the  time) .  At  every  variation  with 
time  of  the  density  of  circulation,  free  vortices  of  the  density  e(x,t) 
will  separate  from  the  lifting  vortices  and  will  drift  away  with  the  air 
flow.   On  the  basis  of  the  theorem  about  the  conservation  of  circulation, 
or  by  integration  of  Euler's  differential  equation  once  along  the  pressure 
side  and  once  along  the  suction  side  of  the  wing  (compare  the  unabbreviated 

"Das  ebene  Problem  des  schlagenden  Flugels"  Zeitschrift  fur  augewandte 
Mathematik  \md  Mechanik,  Band  \,    \S2}\,   pp.  277-292. 

Abstract  from  the  Gottinger  dissertation  of  the  same  title.  Avail- 
able in  the  University  library  at  Gottingen  and  the  state  library  at 
Berlin.  Referent:   Professor  Prandtl. 

^This  periodical  (Z.f .a.M.M.) ,  vol.  3,  1923,  pp.  290-297. 
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report) ,    one  finds  easily  that  the  vortices  must  satisfy  at  every  point 
and  at  every  time  the  following  continuity  eq.uation  of  the  vortex  density-^: 


^  +  M  +  V  ^  =  0  (1) 

St   bt  Sx 


Generally  one  will  make  suitable  assimptions  regarding  7(x,t)   so  that 
e  may  be  found  from  the  equation.   If  by/dt     is  designated  by  7(x,t), 
its  integral  is 


e(x,t)  =  -  ^  /       7(1, t  +'i^)d|  (2) 

""^^Cx-vt)   ^       ^     ' 


Where  9  is  an  arbitrary  function  to  be  determined  by  bovmdary  condi- 
tions. After  e  has  thus  been  fo\ind,  there  results  the  induced  vertical 
velocity  in  first  approximation  as  principal  value  of  the  integral 


w(x,t)  =  J-  P  2_L1  d|  (3) 

2n  J_   X  -  F 


Similarly  to  the  former  case,  w  now  has  to  be  put  in  relation  with  the 
moved  wing  contour.   If  the  motion  of  the  wing  is  indicated  by  the  value 
of  the  ordinate  at  the  point  x,  thus  by  y  =  y(x,t)   (with  this  formu- 
lation one  could  also  include  a  moved  and  simultaneously  deforming  wing) , 
the  kinematic  boundary  condition  reads 

^  +  V  ^  =  w(x,t)  iS) 

dt     dx 

■^Details  on  the  derivation  of  this  equation  may  be  found  in  a  lecture 
of  Prandtl  "On  the  Formation  of  Vortices  in  the  Ideal  Fluid"  at  the  hydro- 
aerodynamic  conference  at  Innsbruck  1922.   (Lectures  concerning  the  field 
of  hydro  and  aerodynamics,  edited  by  Th.  v.  Karman  and  T.  Levi-Civite, 
Berlin  I92U)  and  in  the  original  report  quoted. 
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This  is  the  equation  from  which  contour  and  motion  of  the  wing  are 


obtained  if 
above  said 


y ,   and  therewith  w,  is  prescribed.   In  analogy  to  the 


i|r(x-vt) 


wU,t  + 


<i| 


(5) 


If  one  takes  into  consideration  that  the  air  velocity  relative  to  the 

supporting  vortex  line  has  the  components  v  and  w  -  r^,    there  result 

ot 

as  before  the  air  forces  for  the  unit  length  of  the  wing 


K 


7  dx 


(6a) 


M  =  pv 


U 


^^1 

1 
'2 


7X  dx 


(6b) 


w  =  p 


w  - 


m 
^t 


dx 


a2 
p.- 


pjr 


(7) 


The  last  term  in  W  represents  as  in  my  first  report  the  suction 
force  at  the  leading  edge  of  the  wing,  now  of  course  with  a  coefficient 
variable  with  time  of  the  first  fundamental  function. 

2.  Periodical  wing  motion.-  From  these  general  expressions,  I  now 
pass  to  treating  the  special  problem  of  the  flapping  wing.   It  is  to  be 
expected  that  7   itself  in  this  case  will  be  periodical,  except  for 
a  constant  mean  value.   Thus,  I  determine  first  from  the  expression 


7  =  7o(x)ei^^  =  7oe^"^^  =  7oe"^'^ 


(8) 
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the  free  vortices  e;  for  7  which  alone  enters  into  e  is  independent 
of  the  constant  mean  value  7  so  that  this  expression  is  sufficiently 
general.  As  custoraa.ry  in  oscillation  theory,  one  tacitly  deals  with 
the  imaginary  constituent  of  7.   v  is  the  circular  frequency  of  the 

V  If 

oscillations,  cd  =  —  a  dimensionless  quantity  which  I  call  reduced 

frequency"  (v  has  the  dimensions  of  an  angular  velocity  since  it  equals 
the  flight  velocity,  measured  with  half  the  wing  chord  as  unit  length) . 
If  one  designates  the  distance  travelled  by  the  wing  after  a  full  period 
as  the  "wave  length"  X,    one  has 


V     03 


The  further  results  are  obtained  by  series  developments  with  respect  to 
the  reduced  frequency  cd  which  have  the  form 


E Cjun'^"'( log  (a)^  m  ^  n 


These  series  converge  satisfactorily  only  for  small  cd  up  to  about 
CD  =  0.1  corresponding  to  a  wave  length  of  30  ving  chords  or  more.  Thus 
I  presuppose  for  the  present  report  slow  quasi-steady  oscillations  of 
the  wing;  this  does  not  imply,  of  course,  that  the  actual  frequency  can- 
not assume  high  values,  too,  if  only  the  wing  chord  is  sufficiently  small. 
If  I  calculate,  finally,  with  oj'  =  io),  the  formulas  obtained,  at  least 
outwardly,  real  coefficients  which  offer  many  advantages  for  the 
calculation. 

I  now  assume  that  parallel  flow  free  from  vortices  prevails  ahead 
of  the  wing.  Then  I  obtain  with  x  =  x  -  l/2  for  e  and  w,  under  the 
presupposition  that  the  "switching-on  process"  of  the  wing  motion  has 
already  run  its  course  (the  motion  thus  has  become  steady),  the  fol- 
lowing expressions: 

e  =  0;   X  <  -1 


e  =  e^  =  a)"e 


^a)'(vt-x)  ("  e'^'^7o(|)d|       -1  ^  x  ^  1 
e  =  ep  =  oD'e'^'(^-x)  j    e'^7o(T)€       x  >  1 


>     (9) 
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J-l   X  -  I        J-1    X  -  I        Ji 


^^: —  <i|    (10) 

X  -  I 


As  in  ray  first  report,  I  determine  7q{>^)    from  the  three  first 
fundamental  functions^ 


''0  =  ^''Oa  ^   ^^b  ^   ^''Oc 


Where  a,  h,      c  now  are  complex  nimibers.   One  has  then  also 


w  =  w^e^        w^  =  aw.  +  hw.,  +  cw^ 
0  0     Oa     Ob     Oc 


The  integrals  (9)  and  (lO)  cannot  be  evaluated  by  elementary  functions. 
The  last  integral  (lO)  in  particular  leads  to  an  integral  logarithm. 
The  Euler  constant  C  which  appears  as  a  consequence  always  occurs  in 
combination  with  log  2  so  that  I  shall  introduce  as  transcendent 
number 


Z  =  log  2  -  C  =  0.11595  •  •  • 

especially  for  the  present  report.   If  one  takes  this  fact  into  con- 
sideration one  obtains  the  w  in  the  form 

^Oa  "=  ^0  ■*"  °^l"^'  "^  °''2'^'  -'■°S  ^'   ■*"  '^■i'^'^   +  a^oo'^log  cjd'  +  agoj'^  +  a£ao'3log  ui' 

(11) 

One  obtains 


^a  =\ 


^  -  X     ^   -  Ji  _  ^2 


1  +  X     'Ob 


=  \|l  -  x2     7q^  =  x^l  -  x2 
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and  corresponding  expressions  with  p^  and  7^^  are  valid  for  Wq^ 

and  Wq  .   The  coefficients  aj^,   Pj^,  y^     are  polynomials  in  x  the 

coefficients  of  which  are  linear  and  rational  in  Z  (regarding  their 
values  compare  the  original  report) .   Besides^  the  appearance  of  log  co' 
shows  that  a  treatment  of  the  problem  with  a  lifting  line  instead  of 
the  supporting  surface  must  fail  since  to  the  transition  to  the  lifting 
line  there  corresponds  a  transition  to  u)'  =  0  (v  increases  arbitrarily 
for  decreasing  wing  chord) .  For  this  reason,  the  theory  of  the  flapping 
wing  here  presented  is  the  simplest  possible  one. 

Bypassing  the  integral  for  y,  we  are  now  concerned  with  determining 
the  wing  motion  in  an  appropriate  manner  and  bringing  it  into  accord  with 
the  expression  for  7.   Since  the  methods  of  the  lifting  vortex  sheet  are 
linear  and  since  a  deformation  oscillation  will  not  be  taken  into  con- 
sideration, shape  and  mean  angle  of  attack  of  the  wing  may  be  disregarded: 
It  is  then  sufficient  to  calculate  the  oscillations  of  a  plane  wing  with 
the  mean  angle  of  attack  0  and  to  superimpose  these  oscillations  as 
small  fluctuations  with  time  of  the  profile  chord  linearly  on  the  arbi- 
trarily prescribed  profile.   For  the  amplitudes  which  are  assumed  to  be 
small,  one  may  disregard  the  fluctuation  of  the  abscissa  of  a  wing  point 
as  small  of  higher  order,  and  has  then  as  the  most  general  motion  still 
possible  (compare  fig.  l) 

y  =  p(t)  +  x9(t)  (12) 

p  and  cp  are  to  be  developed  as  Fourier  series  with  respect  to  cdv; 
I  retain  only  their  first  term 

y  =  (A  +  Bx)^'^  (15) 

Higher-harmonic  oscillations  would  again  have  to  be  superimposed  linearly. 
Fluctuations  in  the  flight  direction  (which  could  be  expressed  by  periodica 
fluctuations  of  v)  also  will  be  disregarded. 

I  am  designating  the  quantities  A  and  B  as  complex  stroke  ampli- 
tude and  amplitude  of  rotary  oscillation,  respectively,  since  equation  (13) 
is  formed  by  combination  of  the  special  translatory  and  rotary  oscillations 

A  =  1,   B  =  0   (a)       and       A  =  0,   B  =  1  (p) 

Thus  7q  and  all  quantities  linearly  connected  with  it  will  be  linear 
and  homogeneous  in  A  and  B.   For  simplification  of  the  calculation, 
the  cases  (a)  and  (p)  may  therefore  be  calculated  separately  and  after- 
wards be  combined  linearly,  for  instance 

a  =  Aa^,  +  Ban  etc. 
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The  complex  circulation  coefficients  are  found  from  equation  (h)    in  the 
following  manner:   a,  t ,      c     are  expressed  in  the  form  of  the  series 


2       ? 

a  =  a  +  a  o)'  +  a  od'  log  oo'  +  a  cu'   +  a.co'^log  o)'  + 


2    2         ^       ^ 
'  log  CD'  +  a/'cu'-'  +  aToa'-'log  od'  + 


ac-cjD '  log 


Yoo'-'log 


aoOD'-^log  0)'  +  aQ(JD'3log3oo'  +  . 


b  =  bo  + 


c  =  Cq  +  .  . 


(11^) 


Equation  (k)    has  in  our  case  the  form 


w(x,t)  =  (Aid'   +   Bco'x  +  B)ve^'^ 


=  fawQ^  +bwQ^  +  cwQ^)e' 


GD'vt 


This  is  to  be  valid  identically  in  x  and  t.   If  both  sides  are 
interpreted  as  series  in  co'^(log  od')'^^  their  coefficients  must  be 
identical,  a,     b,   c  therein  are  unknowns.   Each  coefficient  yields 
a  relation  in  x  which  at  first  cannot  be  expected  to  be  identically 
f ulf illable  by  only  three  free  values .   It  was  intended  to  fulfill  the 
condition  at  least  at  three  suitably  selected  locations  X]_,   Xg,  and 

Xo.   Surprisingly,  these  conditions  are  now  shown  to  be  precisely  of  the 
second  and  not  of  a  higher  degree  in  x.  At  least  this  is  the  case  up 
to  the  terms  of  third  order,  and  there  is  no  doubt  that  this  will  remain 
so  for  the  higher  terms  as  well.   The  coefficients  of  second  degree 
in  X  determine,  therefore,  with  their  three  subcoeff icients  each  one 
triplet  of  values,  each  of  the  a,  b,   c.   In  every  new  comparative 
coefficient  there  appears  a  new  value  triplet  of  this  kind  so  that  the 
thirty  imknowns  for  which  a  formulation  has  been  set  up  may  be  foimd 
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uniquely  from  linear  relations  by  successive  evaluation.  The  calculation 
results  in  the  following  values: 

ao  =  2Bv 


ai  =  2vfA  +  ^  B(1  -  2Z)] 


ag  =  2Bv 


[Taz  -  b(z  -  z2)] 


a^  =  2v  l-AZ  -  B 


ai^  =  2v[a  +  B(1  -  2Zy| 


a^  =  2Bv 


ag  =  2v  AZ2  -  Bfi  +  i  Z  -  |  Z^  +  z3   (See  footnote  5.) 
a^  =  2v  p2AZ  +  B(i  -  Z  +  J^zA 
ag  =  2v[a  +  B^l  -  3Z)] 


an  =  2Bv 


bo  =  0 

b-j_  =  Ij-Bv 
b2  =  0 


b,  =  2v(a  -h  I  b) 


3 

Cq  =  0 

c-L  =  0 

C2  =  0 


Co  =  Bv 


b^  =  Ci  =  0,   i  ^  Ij- 


The  dimensions  of  all  quantities  are  affected  by  the  selection 
of  the  wing  chord  as  unit  length.  A  treatise  of  the  author  in  the 
Zeitschrift  flir  Motorluftschiffahrt  on  the  same  subject  has  been  arranged 
so  that  no  objections  are  possible  from  the  viewpoint  of  similarity 
mechanics . 
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I  want  to  point  out  here  briefly  that  approximated  values  for  the 
circulation  coefficients  are  obtained  also,  if  the  effect  of  the  free 
vortices  is  disregarded  and  the  elementary  calculation  made  in  such  a 
manner  as  if  the  momentary  apparent  angle  of  attack  of  each  wing  element 
were  decisive  for  its  circulation,  namely  the  angle  formed  by  the  air 
velocity  v  relative  to  the  moved  wing  element  and  the  direction  of  the 
latter.   Since  in  case  of  rotary  oscillations  (B  j^   O)  every  wing  element 
has  another  vertical  velocity,  the  apparent  angles  of  attack  of  the 
elements  are  all  different,  that  is,  the  wing  assumed  to  be  plane  behaves 
as  if  it  had  an  apparent  (dynamic)  cuin/-ature  which  is  periodical.  The 
simple  calculation  (compare  the  original  report,  third  part-,  beginning 
of  section  II)  yields  accordingly  circulation  contributions  of  the  two 
first  fundamental  functions,  that  is,  only  the  following  terms: 


aQ  =  2Bv 


bQ  =  0 


i^  =  2v(a  +  I  b)      b^  = 


2Bv 


Thus  except  for  higher  terms  and  with  consideration  of  the  order  of  magni- 
tude of  B  (see  below)  a  good  approximation  is  obtained. 

All  the  rest  follows  readily  from  the  circulation  coefficients.   I 
had  introduced  the  quantity  od'   in  order  to  enable  an  easier  calculation 
and  consideration  also  of  complex  <d',  that  is,  damped  or  excited  oscil- 
lations.  It  is  true  that  for  damped  oscillations  the  integrals  (lO)  lose 
their  meaning,  since  the  oscillation  had  been  assumed  to  have  been  going 
on  for  an  infinitely  long  time  so  that  here  infinitely  large  amplitudes 
would  have  had  to  precede.   However,  the  formulas  obtained  may  be  retained 
as  approximations,  provided  a  correction  is  made  for  the  starting  process. 
For  the  following  consideration  of  the  steady  oscillations,  I  continue  the 
calculation  for  reasons  of  physical  clarity  with  the  real  frequency  00. 

Lift  and  moment  are  purely  periodical  quantities,  that  is,  their 
temporal  mean  value  is  0.  Nevertheless,  their  amplitude  and  their  phase 
angle  are  of  interest.  According  to  formula  (6),  there  is 


K  =  pv^jr^Akot^  +  Bkp)e 


ioDvt  _ 


M 


pv2jt(AkQ^  +  Ckyjei'^^ 
pv^jtcorASo^  +  Cky)ei'^>vt 


=  pv^rt^Am^^  +  Bm3)e^'^^^  =  pv^it^AmQ^  +   Cm~)e 

=  pv^jtcjD^Aii^  +  CiSyj  e 


.iojvt 


iODVt 


(15) 


10 
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Therein 


B  =  Coj;   k^  =  oikpj   m^  =  ojnp;   k^^  =  dk^', 


1%  =  'iUTa;   kp  =  kyj  mp  =  HLy 


I  have  introduced  here  the  quantity  C  as  new  amplitude  of  rotary  oscil- 
lation. This  was  done  because  in  case  of  ordinary  flapping  motions 
B  is  of  the  order  of  magnitude  oA  as  a  simple  consideration  shows. 
The  value  of  00  is  assumed  to  be  small;  thus  A  and  C  will  be  of 
the  same  order  of  magnitude  so  that  calculation  with  C  instead  of  B 
will  be  more  convenient  in  practice.  The  third  form  of  the  air  forces 

permits,  for  cases  of  equal  stroke  velocity  Aioie-'-^^,  comparison  of 
these  forces  in  a  simple  manner;  the  roughest  approximation  theory  wovild 
yield  constant  coefficients  k^,  =  2i,  k^  =  2.  The  coefficients  k  and 
m  are  complex  numbers  the  constituents  of  which  are  given  by  the  fol- 
lowing series  developments .   It  is  noteworthy  that  the  series  for  m 
are  finite. 


^a  =  '^0'   -^  i^a"»     %  =  "V'   +  "\j"'     <^  =  '^'^'7 


c^" 


a^'    =  k^'    =  -(1  -   2Z)a)^   -   2a)  log  od  -   ZnZay'  +  2jta)3log  cd  + 


=  ka_"  =  2a)  -  na)^  +  V^  "  ^^J^^  +  ^2a)3log  o) 


>  (16) 


20)^  log  0)  + 
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kg'=ky'=iky'=2-na)+  i^l  +  2Z  -   2Z^  -  iW^ 


2(1  -  2Z)ao2log  CO  -   2co^log^a3  -    (l-  -  2.  +  nZ  -  ^nZ^juy'  + 


jt(l  -  6z)aPlog  0)  +  3i^<^  log  ^  +   • 


"    _  V    "   _  1 


=  k  "  =  i  ky"  =  (3  -  2Z)a)  +  20)  log  00  -  rt(l  -  2Z) 


ay 


2n:cjo  log  cjd  + 


lit.  I 
\h        2 


+  Z  -  I  jt^^  _  z^  +  2z3)co3  + 


(I  jt^  -  1  +  2Z  -  ez^jcD^log  CD  -   (1  -  6z)a)^log 


CJD    - 


2aPiog3cD  + 


>       (17) 


ma,'   =  aria,'    =  "  2  "^^ 


(18) 


mp '    =  nty '   =  j5  m^ '   =  -  ^ 


3  ,,.2 


CD 


nin 


™    "         1   ™    ' 


=    CD 


(19) 


It  suggests  itself  to  represent  k^  and  m^j  by  the  initial  posi- 
tion of  "time  vectors,"  visualized  as  rotating,  in  the  complex  number 
plane  (as  customary  in  alternating-current  techniques)  and  to  combine 
from  them  -  with  the  parameters  A  and  C   (of  which  A  may  be  assumed 
real  without  impairing  the  generality)  -  linearily,  in  the  known  manner, 
the  ajnplitude  coefficients  k  =  Aka,  +  Cky  (and  correspondingly  m) 

according  to  magnitude  and  phase.  The  diagram  (figs.  2  and  3)  shows  the 
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curves  of  the  end  points  of  the  vectors  k  and  m  as  functions  of  the 
parameter  cd.   In  the  representation  of  the  cvtrves  for  k  and  m  which 
would  yield  the  most  accurate  values  for  the  graphical  evalviation,  the 
curves  of  the  various  approximations  have  been  plotted  side  by  side  for 
comparison  of  the  convergence  of  the  series .  Of  particular  interest  is 
the  case  where  the  wing  -  without  being  affected  by  significant  air 
forces  -  glides  over  an  undulating  streamline  course,  clinging  to  it 
as  much  as  possible.  To  this  corresponds  the  parameter  C  =  -iA,  which 
in  fact  yields  small  air  forces  of  second  order  in  o),  namely 

i  k'  =  2^0^  +  .  .  .       i  k"  =  [i  -  2Zja)3  +  .  .  . 

i  m'  =  i  cd2  i  m"  =  a  0)3 

A      2  A      8 

3.  The  induced  drag.-  The  induced  drag  is  no  longer  linear  in  the 
circiilation  so  that  the  complex  method  could  not  be  retained  without 
new  stipulations.   It  offers  no  longer  any  simplifications,  and  it  is 
advisable  to  continue  from  here  on  the  calculation  with  the  imaginary 
constituent  of  all  quantities  in  real  form.  If  the  calculation  is 
carried  out  according  to  equation  (7),  W  assumes  the  form 

W  =  Wq  +  W-L  sin(ajvi:  +  q)j_)  +  W2  sin(2(jmrt  +  cpg) 

Here  W]_  is  different  from  zero  only  when  the  oscillation  is  superimposed 
on  a  constant  angle  of  attack  different  from  zero.  W  is  in  00  small, 
of  second  order.  The  piirely  periodical  terms  are,  therefore,  hardly  sig- 
nificant; however,  the  temporal  mean  value  Wq^  which  is  different  from 
zero  is  important.  I  calculate  only  this  value  and  write  for  it,  for 
reasons  of  simplicity,  again  W.   I  equate 

A  =  A'  +  lA",   B  =  B'  +  iB",   C  =  C  +  iC" 

and  may  assiome,  without  restricting  the  generality.  A"  =  0.  Then  W 
becomes  a  quadratic  form  in  A'  =  A,  B',  B"  (or  A,  C,  C"):  a  simple 
deliberation  shows  that  the  coefficients  of  B'2  and  B"2  (or  C'2 
and  C"2)  are  equal  and  that  the  coefficient  of  B'B"  (or  C'C")  is 
zero.  Thus  W  becomes 


W  =  pvSnJA^w^  +  2AB'w^p'  +  2AB"w^p"  +  (B'2  +  B"2)wpp| 
«|a2w^  +  2AC'w^y'  +  2AC"w^y"  +  (C'2  +  C"^)v^2. 


2 


(20) 
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The     Wjj^     are  again  series   in     cjuP(cu  log  u))™     and  have  up  to  higher  terms 
the  values 


v^  =  -CO' 


i2  +  jta)3   _  /I  jt2   -  Z^joD^  -   2Z6D^log  ao  +  OD^log^o)  +    .    .    , 


V, 


ap 


w. 


CD 


0.7 


i  (3   -  2Z)cjo2   -  i  oi^log  o)  + 
4     -^  2 


i  (3  -   2Z)(xr'  +  i  cD^iQg  (j^  + 
!<.  2 


^ap     -  ^  ^a7 


1  3        2 

2  1^ 


fl  «2   .  1  ,2  ^  1  ,\ 
V8  2  2     / 


3 


i  (1  -  2Z)co3log  00  =  i  a>3log^a)  + 


1   ^..2 


>(21) 


'3P  -  :j  ^77  -  -^  '^  -  ^  -'   ■  ^'"  '8 


(n^  +  1)cd2  +  1  (3Tt2  +  2   -  Uz2  -  4z)co3  + 


00'- 


i    (1    +    2Z)cD3log   CJO    -    i  (JO^log^OD   + 

2        /    B     2 


W  may  be  positive  or  negative,  according  to  the  selection  of  parameters. 
Depending  on  the  type  of  motion,  one  has,  therefore,  to  expect  drag  or 
thrust.   I  postpone  detailed  discussion  until  after  calculation  of  the 
power  requirement  and  the  efficiency.  The  case  of  gliding  over  an 
imdulating  flow  mentioned  above,  corresponding  to  C '  =  0,  C"  =  -A, 
gives 


W  =  A^pv^jt  -cD^log  o)  -  (i  -  zU^  >  0 


Thus  the  selection  of  parameters  made  does  not  yet  correspond  exactly  to 
the  case  W  =  0;  a  small  correction  would  have  to  be  provided  for  this 
case. 
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k.   Work  done  at  the  ving.-  In  the  free  vortices  hehlnd  the  wing, 
energy  is  contained  which  must  be  produced  by  mechanical  work  on  the 
airplane.  This  may  be  done  in  two  ways.  The  flapping  motion  may  resvilt, 
as  mentioned  before,  in  positive  or  negative  thrust.  In  the  case  of  neg- 
ative thrust  a  propeller  which  overcomes  this  and  all  other  resistances 
to  flight  is  required  for  maintenance  of  equilibrium  of  motion.   In  the 
case  of  positive  thrust  a  propeller  is  needed  only  until  the  thrust  due 
to  flapping  exceeds  the  resistances  to  flight,  whether  the  flight  be 
uniform  or  accelerated.  As  to  the  work  performed  at  the  wing  itself,  the 
wing  motion  consumes,  of  couse,  energy  if  thrust  exists;  the  ratio  of 
thrust  power  -L^.  =  -Wv  and  the  total  mechanical  power  L^  to  be  applied 
to  the  wing  may  then  be  denoted  as  aerodynamic  efficiency  of  the  flapping 
wing.   In  case  of  drag,  two  more  possibilities  exist.  First,  the  flapping 
motion  may  require  additional  work .   The  efficiency  defined  above  then 
becomes  negative  and  arbitrarily  large  when  the  wing  power  L^  decreases 
more  and  more.  Second,  the  case  may  occur  that  Lf  becomes  negative, 
that  is,  the  wing  then  is  supplied  with  energy  from  the  air  (indirectly 
by  the  propeller)  and  may  vise   that  energy  for  surmounting  the  resistances 
in  the  oscillation  mechanism,  or  may  store  it  in  the  oscillation  itself, 
that  is,  increase  its  amplitudes.  Aside  from  this  "incremental  power" 
the  propeller  must,  of  course,  in  this  case  yield  additionally  the  energy 
of  the  free  vortices  so  that  one  may  define  the  quotient  of  -Lf  and 
the  total  propeller  power  as  efficiency  referred  to  the  power  absorption 
of  the  wing.  This  is  then  exactly  the  reciprocal  value  of  the  efficiency 
defined  above  which  in  this  case,  as  the  quotient  of  two  negative  numbers 
becomes  positive  but  larger  than  one,  thus  loses  its  physical  meaning. 
Lf     is  divided  into  two  parts .  One  has 

-^  =  p  +  9X  =  ia3v(A  +  Bx)e^"^ 
ot 

Therewith  L-^  =  Kp  becomes  the  "flapping  power"  and  L^,  =  Mcp  the  power 

opposed  to  the  rotary  oscillation;  thus  the  wing  power  is  Lf  =  L^  +  Lr* 
The  power  opposed  to  the  drag  is  denoted  by  I^  =  Wv.   In  the  air  there 
then  remains  in  all  L=Lf+Lyr  =  L^+Lj.  +  Ly.  I  indicate  of  all  Lj^ 

again  only  the  temporal  mean  values  for  which  I  obtain  qioadratic  forms 
of  the  same  type  as  for  W 

Li  =  pv^nJA^Zi^  +  SAB'Zi^p,  +  2AB"Zi^pM  +  (B'2  +  B"2)  lippj 

y  (22) 

=  pv^nh^l^^   +  SAC'Zi^^.  +  2AC"Zi^^,  +  (C'2  +  C"2)  Z^  1 
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L]-  has  finite  series  and  is,  in  general,  small  compared  to  L+.  ike 
coefficients  are  individually 


^wik  ~  '"'ik 


(25) 


^toa  =  "^^  "  f  ^^  "^  (tT  "  ^^1^^  ■*"  2Za)^log  oj  -  cjo^log^cu  + 


Ha3'    =  I  Ha/'    =  J  (3   -  2Z)a32  ^  |  cD^log  co 


i   (1   -    2Z)od3    -  ^  OO^log  ^3  + 
4  2 


•,  1   ■,  1  Jt      2        /it'^        1  „        1  ,,2        1\    3 

tap  ta7         2  U  V8        2  2  8/ 


i  (1  -   2Z)cL3log  0)  -  i  CD^iogS^^  + 


Hpp  -  ^  H77  -  ° 


>       (2i|) 


^roa  =  0 


Wp'   -  ~  Va7'   -  ° 


1  1     ^ 

^v~^q"    =  —    ira-y"    =  rr  cjD-^ 


rap' 


CO 


''ra/ 


7  -    1     7  _  1     2 

^rpp  -  ^   ir77  ~  2  ^ 


)- 


(25) 
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^foa  =  0)2  .  I  0.3  + 


{i  -  ^^) 


4  i4-  k         9 

o)     +  SZfOO  log  CD  -  o)  log  O)  + 


O) 


^fap-   =  ^  ^far'    =  ^  (3  -  2Z)a)2  +  i  a)2iog 


21  (1  _   2Z)ol3   -  i  a)3log  OD  + 
4  2 


IfaB"   =  -   ifa-x"   =icD-«a)2+(4  +  -Z-i  22^5 
lap         CO     far  2  k  V8        2  2      / 


i  (1  -   2Z)a)3log  CO  -i  ooSlog^o)  + 
2  2 


1  1     1  12 

or 


>        (26) 


3t        o  n''        k 

cux       2  2 


1  ■^ 

Z„o  I   =  —  Z_-y  I   =  —  oi-^  + 


1  ,  «     2       n^     3    ^ 

II     =    _   CD       -    00-^     + 


^ap"  -  -  ^ar         2 


^PP  =  ;^  Vr  =  f  - 


Hi  cd2  W I  «3  ^  i 

2  l8  4 


2L  Z'^ 
2 


21  zjcD^  + 


-    (1    +    2Z)£D3log   CD    -    -  Co3log2(D    + 


(27) 
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L,  as  the  energy  of  the  vortex  trail,  can  of  course  never  become  nega- 
tive and  must  therefore  "be  a  positive  q.\iadratic  form  -  definitive  in  the 
amplitudes .  The  fact  that  L  in  the  form  here  noted  is  capable  also  of 
small  negative  values  is  caused  by  the  neglect  of  higher  terms  in  the 
series  development.  What  is  obtained  in  this  case,  is  therefore  only  the 

error,  accidentally  negative,  of  the  almost  vanishing  vortex  power. 

p 
Since,  when  C   is  used,  all  coefficients  contain  the  factor  o)  ,  the 

latter  has  been  cancelled  in  the  graphical  plotting  which  is  in  agree- 
ment with  the  presentation  of  the  power  for  constant  flapping  velocity 
(fig.  k)  .      The  c\arves  show  that  the  essential  terms  always  stem  from  the 
stroke  amplitude  A,  possibly  in  combination  of  the  latter  with  the 
amplitude  of  rotary  oscillation,  and  that  the  corresponding  coefficients 
increase  somewhat  more  slowly  than  o)^.   By  rotary  oscillation  I  meant 
above  an  oscillation  about  the  z-axis.  More  generally,  every  oscillation 
where  the  ratio  of  A  and  C  is  real  is  a  rotary  oscillation  about  the 
fixed  axis  with  the  abscissa  a  where  A  then  is  A  =  -aCoj.  It  is 
shown  that  for  not  too  large  values  of  am,  that  is,  for  axes  which  do 
not  lie  at  too  great  a  distance,  W  and  L  are  always  positive;  that 
is,  it  is  not  possible  to  obtain  thrust  or  power  absorption  by  rotary 
oscillations  about  fixed  axes.  Production  of  thrust  always  requires  a 
stroke  amplitude  different  from  zero,  power  absorption  req_uires  addi- 
tionally a  rotary  oscillation  lagging  by  about  90°-  Pure  stroke  oscil- 
lation without  rotation  also  produces  thrust  which  in  roughest  approxi- 

2   2   2 
mation  results  as  -W  =  A  pv  nco  ,  similar  to  the  so-called  Knoller-Betz 

effect  (if  one  calculates  with  the  y-axis  as  the  "polar"  in  the  plane 

problem) . 

One  obtains  good  insight  into  the  variation  of  drag  and  power  if 
one  varies,  for  fixed  absolute  value  of  the  amplitude  ratio  C/A  =  c, 
only  the  phase  angle  cp  between  the  two  oscillation  components  where 
one  then  has  to  put 


C  =  Ac'  =  a|c|  cos  cp,      C"  =  Ac"  =  A[c|  sin  cp 


W  is  shown  to  become  a  minimum,  the  thrust  thus  a  maximum,  when  the 
rotation  leads  by  somewhat  more  than  90°,  namely  by  cp  =  arc  tan  — —. 

It  is  plausible  physically,  too,  that  the  thrust  will  assume  large  values 
precisely  then  when  the  phase  is  shifted  by  about  l80°  compared  to  the 
phase  which  is  present  for  gliding  free  from  air  forces  over  the  wave 
course.   For  every   |c|   there  exists  an  00  and  vice  versa  for  which 
the  thrust  maximum  is  absolute.  One  then  has 
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^77 


W   ti 
c  (O))  = 


^min  =  A^pv^rtWj^in 


^77 


=  _ia.-  J-a.2_  J_(l+  3„2)^3  _ 


2jt2      2rt3 


L  becomes  a  minimum  and  disappears  for  suitable  co  except  for  terms 

of  the  fifth  order  when  the  rotation  is  lagging  by  somewhat  more  than  90°. 

Finally,  the  wing  power  L|.  becomes  a  minimum  -  thus,  the  absorbed  power 

a  maximum  -  when  the  rotation  lags  by  somewhat  less  than  90°;  namely  for 

^far" 
cp  =  arc  tan  - — '—.     This  maximum  too  becomes  absolute  when  between  c 

^fa7' 
and  CD  the  following  condition  is  satisfied 


^f77 


c'  (cc)  =  - 

ifyy 


^fmin  =  ^^P^^-^l-fxaln 

Va7'   ^  ^far"^        rt 

ifmin  =  ^foux =  -1  +  2  CO  + 

'f77 

Altogether,  L  becomes  negative  only  when  |c|  >  1. 
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The  efficiency 


,      r        J  T)-]^  is  valid  for  Lf  >  0  and 

T^,  =  —  =  -  —  =  1  -  —  where 

^2    Lf      Lf        ^2     ^°^    ^f  "^  °- 


>      (28) 


is  the  quotient  of  two  q_\iadratic  forms  and  capable  of  a  great  many  values . 
Generally,  it  "becomes  negative  for  small  or  completely  vanishing  stroke 
amplitudes,  and  arbitrarily  large  with  uj — ^0.  The  same  is  valid  for 
rotations  about  a  fixed  not  too  remote  axis.   If,  however,  k  f  Q     and 
also  lim  A  j^  0,   t]  approaches  for  o)  — >•  0  to  the  limit  1. 
(0=0 

By  a  rotary  component  (c"  =  O)   of  equal  or  opposite  phase  tj  is 
always  deteriorated  compared  to  c  =  Oj  the  same  is  true  for  c '  =  0, 
c"  >  0.  In  contrast,  the  efficiency  is  improved  for  c '  =  0, 
-1  ^  c"  ^  0,  as  can  be  seen  from  the  diagram  (fig.  5).  For  c'  =  0, 
c"  =  -1,  T]  becomes  identically  1,  and  for  c'  =0,  c"  <  -1  there 
results  power  absorption.  The  representation  for  fixed   |c|   in  depend- 
ence on  cp  as  a  discontinuous  single-wave-harmonic  f\inction  is  very 
graphical  for  the  efficiency  as  well.  Figure  6  shows  clearly  at  what 
phase  angles  the  transition  from  power  absorption  to  power  production 
takes  place.  The  most  important  ones  among  the  coefficients  found  from 
the  series  developments  have  been  compiled  in  the  numerical  table 
(table  I) . 

5.  Application  to  the  flutter  of  elastically  supported  wings.-  The 
derived  laws  could  be  practically  applied  in  the  investigation  of  a 
phenomenon  our  pilots  observed  in  the  last  war.  In  the  so-called  sesqui- 
planes,  the  lower  wing  was  fastened  to  one  single  spar  only,  thus  was 
only  slightly  elastic  against  small  deflections  and  rotations.   In  case 
of  increased  flight  velocity,  for  instance  in  steep  dives,  there  occurred 
sometimes  vigorous  flutter  of  the  lower  wing  tips  which  underwent  obviously 
unstable  oscillations  in  the  increased  air  flow.  Of  course,  such  unstable 
oscillations  are  possible  only  if  energy  is  supplied  to  the  oscillating 
system,  and  tMs  occurs,  according  to  ray  investigations,  only  when  the 
vector  of  the  ajtiplitude  of  rotation  lags  by  about  90°  and  when  the  ampli- 
tude of  rotation  itself  is  sufficiently  large  i\c\      must  be  >  |A|j.   Let 
us  visualize  again  the  gliding  -  almost  free  from  air  forces  -  of  the  wing 
over  an  undulating  flow  course  where  the  airspeed  (relative  to  the  wing 
elements)  has  no  vertical  component.   If  the  amplitiode  of  rotation  is 
smaller  than  corresponds  to  this  case,  the  motion  is  damped  by  the 
counteracting  air  force.  If  the  amplitude  of  rotation  is,  on  the  con- 
trary, larger  than  in  the  case  above  and  the  wing  therefore  scoops  more 
deeply  into  the  air,  the  air  force  always  acts  in  the  direction  of  the 
motion,  and  the  motion  is  excited. 


20  MCA  TM  1361^ 

I  consider  a  wing  supported  on  spars,  elastic  with  respect  to  trans- 
lation in  the  y-direction  and  with  respect  to  twist.   In  order  to  be  able 
to  go  on  from  my  formulas  used  so  far,  I  introduce  the  directional  forces 
per  unit  length  in  z-direction;  I  calculate  therefore  as  if  these  forces 
were  distributed  continuously  over  the  length  of  the  wing.   This  assump- 
tion does  not  lead  to  any  contradictions  if  the  wing  in  itself,  aside 
from  its  support,  is  sufficiently  stiff.   Since  I  disregard  deflections 
in  the  flight  direction,  I  can  show  that  a  wing  supported  on  spars  always 
has  only  three  essential  elasticity  parameters,  corresponding  to  the 
three  constants  of  the  work  of  deformation  quadratic  in  p  and  cp . 
(Compare  fig.  1  and  the  original  report.)   With  respect  to  its  elastic 
properties,  this  wing  may  therefore  always  be  replaced  by  a  wing  which 
is  supported  only  on  one  spar  (the  "elastic  axis")  with  the  abscissa  a, 
elastic  with  respect  to  translation  by  means  of  the  directional  force  c, 
and  with  respect  to  rotation  by  means  of  the  directional  moment  7,  as 
schematically  indicated  in  figure  1.  The  resultant  of  the  elastic  forces 
and  its  moment  at  the  origin  then  are 

K  =  -c(p  +  acp) 

M  =  -cap  -  (ca^  +  7)<P 

If  s  is  the  abscissa  of  the  center-of -gravity  axis  (point  S, 
fig.  1),  the  momentum  and  its  moment  at  the  origin  have  the  values  (with 
the  mass  m  per  unit  length  and  the  corresponding  moment  of  inertia 
■iS  =  mr^  for  the  center -of -gravity  axis) 

G  =  m(p  +  scp) 

U  =  rasp  +  (ms^  +  i3)cp  =  msp  +  m(s^  +  r^)cp 
The  eq.uations  of  motion  of  the  wing  then  read 

G  =  K      U  =  M 
For  K  and  M  the  air  forces  have  to  be  added.   If  I  put  as  before 

p  =  Ae"^'^      cp  =  Be"^'^ 
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and  take  into  consideration  that  the  air  forces  have  on  the  wing  the 
opposite  effect  from  the  one  the  forces  of  the  wing  calculated  above 
have  on  the  air,  K  and  M  become 

K  =  -e'^'^  jcA  +  caB  +  pv^itCAl^  +  Bkp)r 

M  =  -e"^'^  -^caA  +  (ca^  +  ■y)B  +  pv^it{Ar[a  +  ^^)\ 
For  abbreviation,  I  introduce  the  following  designations 


c      2 


mv 


7    2 
c    ^ 


Therewith  there  result  finally  as  the  equations  of  oscillation 


A(od'2  +  0^2  ^  p^^j  ^  B(su)'^  +  ax^2  +  p^j^.^)  ^  q 
A(sa)'2  +  aacj-,^  +  PQm^)    +  bUt^  +   s2)a)'2  +  (a^  +  ^2)^:^2  +  PQmp^  =  0   (29) 


Equating  the  determinant  to  zero  yields  the  eqimtion  for  the  fre- 
quencies OD  =  -im' .  If  I  first  disregard  the  forces  (for  this  purpose 
I  put  P-.  =  0) ,  there  result  two  main  oscillations  as  a  consequence  of 

the  coupling  of  the  oscillation  of  the  center  of  gravity  with  the  fre- 
quency a:^-,  and  the  oscillation  about  the  center  of  gravity  with  the 

^  /,^  i.r-i+v>      p,_^o    _o\2j.^2j.p,2 

r 


frequency     -  oc^.     With     5   =   (s   -  a)^  +  r2  +  q2     there   i£ 


03.    .2  .   ^^    ^.2  ^  IU_^g   ^  \|52  ..  I.q2r2) 


,2  .'^ 


^1.2 


1.2 


2r' 


A 
B 


2ar' 


1^6   ±  ^6^   -  i^q2r2j 


&   ±  \S^   -  i4-q^r2   -   2r2 


=  -a 


1.2 


y^  =  B^{x  -  ak)ei'%^j      k    =  1.2 


(30) 
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The  main  oscillations  therefore  are  rotary  oscillations  about  fixed 
axes  at  the  distances  a^^ag.  The  quantity  s  -  a  forms  a  measure  for 

the  coupling.   In  general,  there  develop  beats  from  both  main  oscilla- 
tions.  If  I  write  those  in  the  form 


y  = 


-a-]_B-]_e 


ievt 


agBge-i^^  +  X 


|B,el 


ievt  ^  B2e-"te^ 


9} 


,ic0\rt 


tD2_  +  CDg  =  2a) 

a>]_  -  (02  =  26  «  0) 


The  motion  may  be  interpreted  as  an  ordinary  oscillation  with  an  ampli- 
tude ratio  slowly  variable  as  to  magnitude  and  phase.   Therein  there 
appears  of  course,  periodically  recurrent,  the  phase  angle  which  corre- 
sponds in  the  air  flow  to  the  power  absorption.   If  the  air  forces  are 
to  counteract  the  change  of  this  phase  angle,  corresponding  to  con- 
tinuous supply  of  energy,  it  can  be  shown  that  the  case  of  slight  coupling 
(s  -  a  small)  for  balanced  or  almost  balanced  frequencies  of  the  uncoupled 
system  is  the  best  presupposition  for  this  phenomenon.   I  am  anticipating 
from  the  results  of  the  following  calculation  with  consideration  of  the 
air  forces  that  without  coupling  (s  =  a)  no  unstable  oscillation  at  all 
woiold  be  possible. 

The  air-force  coefficients  occurring  in  the  oscillation  determinant 
are  themselves  functions  of  the  frequency  o)'  which  cannot  be  indicated 
by  simple  analj'-tical  expressions  so  that  the  roots  of  the  oscillation 
determinant  cannot  be  obtained  in  a  simple  manner.   However,  since  their 
existence  is  secured  by  the  physical  meaning  of  the  problem,  it  is  per- 
missible to  introduce  into  the  eqnoation  instead  of  the  coefficients  k 
and  m,  the  first  terms  of  their  series  developments,  all  the  more  so 
since  the  occurring  factor  pQ  generally  is  a  small  number;  for  it  is 

sensible  to  break  off  the  series  for  cos  x  in  order  to  find  from  the 
polynomial  obtained  for  instance  approximately  the  first  zero  of  this 
function  whereas  the  same  method  is  meaningless  for  e-^,  since  no  zeros 
of  this  function  exist  in  the  finite  domain. 

The  oscillation  determinant  obtains  the  following  form 


r^cD'^  +  Scjuq^o)'^  +  q^a:^3^  +  PqCd'^  (r^  +  s^)k^   -  sko  -  smj^^  +  mT]  + 
Po'^O^il^^  ■*■  ^^^^a,   -  ^p  -  a^  +  nip]  +  Po^(^"ip  -  ^p°t(,)  =  °       (31) 
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The  roiighest  approximation  is  the  result  of  the  calculation  which 
bases  the  determination  of  the  lift  and  moment  coefficients  -  in  the 
manner  of  the  theory  of  the  Knoller-Betz  effect  -  on  the  momentary 
apparent  angle  of  attack  and  the  apparent  (dynamic)  curvature  of  the 
wing.   Even  with  this  procedure,  there  result  in  certain  cases  complex 
roots  gd'  with  positive  real  parts  which  resiilt  in  an  "increment"  of 
the  oscillations  and  thus  correspond  to  dynamic  instability.   I  shall 
not  here  discuss  this  approximation  more  closely.  Also,  I  shall  only 
briefly  mention  the  case  of  the  greatest  instability  since  this  case 
is  trivial.  It  occurs  if  the  elastic  axis  lies  so  far  to  the  rear  that 
at  the  slightest  displacement  of  the  wing  from  equilibrium  position  the 
air  flow  simply  causes  the  apparatus  to  tip  over  aperiodically  toward 
the  rear.  This  always  occurs  as  soon  as  a  is  positive  and  the  air 

velocity  sufficiently  large,  namely  for  a  >  or  v^  >   

2pQ  2pjta 

for  a  >  0. 

If  I  retain  of  k  and  m  all  terms  up  the  second  order  in  go', 
the  period-equation  of  the  wing  becomes 


f(a)')  =  aQU)'^  +  aQ'oj'^log  co'  +  a^'ui'^log^cD'   +  bQOo'^  +  bQ-co'Sipg  cu'  + 

CqU)'^  +  CQ'CD'^lOg  U)'  +  CQ"a)'^l0g^CD'  +  dQGD'  + 

dQ'oo'  log  CD'  +  Cq  =  0  (32) 


With  the  coefficients 


ao  =  r^  +  p^6  + 


PO^   -^  Po'(i  +  i  Z  -   Z2)   >  0 


bo  =  Pq^  +  Pq^(^  -  z)    >  0 


Cq    =    CDQ^^g     +    p^v)      -     2PqS     +    Pq^ 


dQ  =  cDq^PqX  >   0 


2k 
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-  ^  2  /  2^  2 
'0  0    V-     0 


en   =  ^n    (l^n      "   SpQa 


^0     =  PQ 


2(r2  +  s2)    -   2s(l  -   2Z)    -   p^^i  -   2z\ 


^0"   =  "PQ^^^   ■"  PO^ 


^0'    =   -Po(2^   -   Pq) 


Cq'    =   2p(^2p   +  q^2   _  g^^   _   2z)J 


Cq     =  do'    =   -2poao^' 


5   =   (s-a)2  +  r2  +  q^2>o 


^   =  1  +  2(r2  +  s2)    _   s(3   -   2Z)   >  0 


X  =  1  +  2(a2   +  q2)    _  a(3   -   2Z)   >  0 


e   =   (r2  +  s2)(l  -   2Z)    -   s(l  -   2Z  +  2z2)    +  |  >  0 

8 


V  =  (q2  +  a2)(l  -   2Z)    -  a(l  -   2Z  +  2z2)    +  1  >  0 

8 


z  =  0.11593 
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6.  Numerical  evaluation.-  The  equation  can  be  solved  only  approxi- 
mately, of  course.  For  this  purpose,  I  first  omitted  the  logarithmic 
terms  and  determined  the  roots  c5'   of  the  algebraic  epilation 


(cjd')  =  agO)'^  +  bg^'3  +  CqCjo'^  +  dQOj'  +  Cq  =  0         (33) 


g 


I  regard  this  equation  as  an  approximation,  equate  cd'  =  co'  +  X,  and 
can  now  develop  the  logarithmic  terms  retaining  the  linear  terms  in  X. 
Thus  I  obtain 


T..  R 


X  =  -od'  log  o)'  I  (3U) 


R  =  a'oi'^  +   a  "cD'^log  00'  +  b  'od'^  +  c  'cu'  +  d  '(l  +  od'  log  cd') 


S  =  imoCu'^+  aQ'a)'5(i  +  U  log  u)')    +  2aQ"aD'3log  a)'(l  +  2  log  o)')    + 

3bQa)'^  +  bQ'aJ'2(l  +  3  log  cB)    +  2cQa3'   +  Cq'od'(1  +  2  log  co')    + 

dg  +  dQ'(l  +  2od'    log  d)')(l  +  log  do') 

The  procedure  may  be  continued  and  yields  the  fiirther  approximation 

R(a3')aL)'    log  o)'   +  g(cD') 


X'   = 


S(a)') 


Finally,  there  follows  the  complex  amplitude  ratio  B:A  =  b  from 
one  of  the  equations  (29) . 
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Finally  one  now  has  to  learn  the  conditions  for  which  the  equation 
for  cjd'   has  complex  roots  with  positive  real  part,  corresponding  to 
excited  oscillations  or  critical  support  of  the  wing.   For  the  limiting 
case  of  dynamic  indifference,  I  formulate  the  roots  of  the  eqiiation  (32) 
as  p\arely  imaginary  and  obtain,  hy  setting  the  real  and  the  imaginary 
constituent  of  the  equation  equal  to  zero,  the  conditions 


aQ  -  —  Slqjci:)     +   a^'cD  log  ao 


+  a^"a)^loe^a)  +  b^ '  ^  o^^ 


0  2 


/  Jt^        „1    2  ,2-,  J    I    /«  -,      2   ^ 

Icq  -  —  c^   IcD     -   Cq  (D  log  (D  -  d_'a3l—  +  (JO  log  00  I 


a_'   -^  (J^     +  aQ"ncjD-^log  03  -  bQ(ja     -  bQ'co  log  co 


Cp, '  T^  CD  +  d^'    log  cd(  1  -  Jtoo)  +  dp,  =  0 


+  e^  =  0 


V     (35) 


'0  2 


0 


From  them  ao  would  have  to  be  eliminated  in  every  special  case 
whereby  a  conditional  equation  for  any  of  the  parameters   s,   a,  r^ , 
1  )     (i^o,      pQ  or  for  the  flight  velocity  v  results.  As  an  approxi- 
mation, it  is  sufficient  to  investigate  the  equation  (33) •   Since  the 
approximation  (3^)  shows  in  the  case  of  instability,  a  small  additional 
damping,  the  criteria  for  instability  derived  from  the  following  equa- 
tions are  necessary  but  not  always  sufficient;  they  are  therefore  valid 
only  with  the  reservation  of  checking  with  the  approximation  (3^) •  When 
they  are  satisfied,  the  wing  may  at  any  rate  be  regarded  as  critically 
supported.   Before  I  set  down  the  most  general  criterion  for  instability, 
I  shall  mention  one  which  is  very  simple  but  suffices  only  for  cases  of 
great  instability.   This  is  the  condition  that  the  coefficient  Cq  becomes 
negative 


v2  > 


prt(2s  -  Pq) 


(36) 
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Written  as  a  function  of  s ,  the  same  condition  reads 


?      2 


2s  (pQ  +  ^) 


Pq2  +  0:^2  /g^2  ^  j,2  ^  ^^2 


Pqv)  <  0 


Fiirther  limits  are  yielded  "by  the  Routh  condition 


^  ^0 


or 


^0  <  ^0  —  +  ^^ 


— —     -r     c-.     

^0  ^'^ 


0^2  (S  +  Pqv)    +  Pq^   _   2PqS    <  |^r2  +  p^e   + 


>(5T) 


H     +     PQ 


e.-) 


In  order  to  satisfy  this  equation,  it  will  above  all  be  required 
that  V  be  sufficiently  large,  that  is,  (Xiq     sirff iciently  small.  Fur- 
thermore,  s  must  be  positive;  therewith  s  -  a  is  positive,  too, 
because,  as  mentioned  before,  a  >  0  would  for  siiff iciently  large  v 
always  result  in  great  instability.   One  can  readily  understand  that 
this  must  be  so.  The  centroidal  axis  of  the  wing  as  axis  of  inertia 
generally  lags  behind  the  elastic  axis  as  line  of  application  of  the 
directional  force.   If  the  centroidal  axis  therefore  lies,  in  flight 
direction,  behind  the  elastic  axis  (s  -  a  >  O) ,  the  wing  has  in  its 
upward  motion,  on  the  average,  a  positive  angle  of  attack;  the  opposite 
is  true  for  the  downward  motion.   Thus  the  air  forces  always  take  effect 
in  the  direction  of  the  motion  and  amplify  the  oscillation;  whereas  in 
the  case  s  -  a  <0  the  opposite,  that  is,  damping  occurs. 

It  can  easily  be  confirmed  that  both  degrees  of  freedom  must  act 
together  for  achievement  of  the  oscillation.  The  calculation  always 
yields  roots  with  negative  real  parts  if  one  of  the  degrees  of  freedom 
is  suppressed  (corresponding  to  ovn  =  00  or  q2  =  00)  .   This  fact  is 

confirmed  by  the  failure  of  tests  undertaken  formerly  in  Gottinger  with 
a  wing  with  only  one  degree  of  freedom. 
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The  power  L  produced  or,  respectively,  absorbed  by  the  wing  is 
according  to  previous  formulas 


L^  =  Prpnv^lf 


(38a) 


^f  =  ^foa  ''■  2^'^fa3'  "^ 


If  I  put  temporarily  co'  =  3  +  io),  the  mean  value  of  the  entire 
wing  energy  for  constant  amplitude  A  is 


—        1         2   2 


Q  =  0^^  +  o)^  +  2b '  iaiar?  +  su?  ]  + 


1 '  Uia^  +  su?j 


(b'2  +  b"2)  \L^  +  <i2)^q2  +    (^2  +  s2)^2^ 


A  =  PiQe^^ ,   and  hence,  in  a  different  form 


_Lf  =  —  =  i  mv^AAQ  =  ^  mvSA^Q 
^   dt   2         2 


(38b) 


If  the  oscillation  calculation  has  been  carried  out,  the  agreement 
of  the  values  l^     found  by  different  methods  offers  therefore  a  control 
for  the  calculation. 

7.  Comparison  with  a  test  result.-  For  confirmation  of  the  theory, 
tests  with  a  light  wooden  wing  of  60  cm  length  and  10  cm  chord  were 
performed  in  the  small  wind  tunnel  of  the  Gottinger  aerodynamic  test 
institute. 
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This  wing  was  suspended  on  an  axis  between  plane  side  walls  by 
means  of  springs.  The  springs  acted  on  cross-shaped  metal  sheets 
which  carried  moreover  sliding  weights  for  variation  of  the  distance 
from  the  center  of  gravity  and  of  the  moment  of  inertia.  The  suspen- 
sion springs  or,  respectively,  their  initial  compression  transferred 
the  directional  force  or,  respectively,  the  directional  moment  corre- 
sponding to  my  formulation  to  the  suspension  axis.  I  shall  not  discuss 
here  fiirther  details  of  the  apparatus;  I  refer  instead  to  the  drawing 
figure  7.  The  purpose  of  the  tests  was  the  determination  of  the  incre- 
ments or  decrements  of  the  oscillations;  for  this  it  was  sufficient  to 
plot  the  oscillation  of  a  point,  for  instance  the  suspension  axis. 
This  was  done  with  the  aid  of  a  registering  drum  kindly  put  at  our 
disposal  by  the  physiological  Institute  of  the  University  Gottingen. 
A  few  of  the  diagrams  of  damped  and  increasing  oscillations  thus 
obtained  are  reproduced  in  figure  8.   From  these  diagrams  the  ratios 

a  =  e    of  amplitudes  succeeding  one  another  were  determined  and  were 
conipared  with  the  values  foimd  by  calculation  with  consideration  of  the 
frictional  damping  produced  by  Oq-j-  and  a^j,  of  the  translatory  or 
rotary  oscillations.  Within  the  considerable  limits  of  error  which 
could  have  been  cut  down  only  by  a  major  expenditure  of  time  and  means, 
the  agreement  may  be  called  satisfactory. 

The  variation  of  the  n\mbers  was  in  excellent  agreement  with  the 
theory.  An  increase  or  damping  of  the  small  oscillations  resulted 
according  to  whether  the  wing  was  overweight  toward  the  rear  or  the 
front;  the  ratio  of  amplitudes  succeeding  one  another  increased  when 
the  overweight  was  increased  or  when  the  air  was  blowing  more  strongly 
against  the  wing.  A  precession  of  the  leading  edge  of  the  wing  was 
clearly  evident  when  the  wing  was  supported  so  as  to  be  unstable .   This 
corresponds  in  the  usual  notation  to  a  lag  of  the  rotary  oscillation 
as  must  be  the  case  for  power  absorption. 

Finally,  it  is  noteworthy  that  for  small  air  velocities  of  about 
5  m  per  second  the  influence  of  the  viscosity  was  shown  by  the  fact  that 
the  lift  was  not  yet  fully  developed  at  small  angles  of  attack.  For  the 
smaller  characteristic  values  of  about  500  m/sec.  mm  a  symmetrical  wing 
has  a  small  range  "dead  angle  of  attack"  where  the  lift  remains  near 
zero  (as  shown  also  by  other  tests  in  Gottingen) .  The  small  oscillations 
were  therefore  still  damped  when  oscillations  with  a  larger  initial 
amplitude  were  already  increasing.   The  limiting  amplitude  lying  between 
these  two  cases  decreased  of  coijrse  more  and  more  with  increasing  v. 


Translated  by  Mary  L.  Mahler 
National  Advisory  Committee 
for  Aeronautics 
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TABLE  I 
NUMERICAL  TABLE  OF  CALCULATED  COEFFICrENTS 


0)  = 

0.02 

0.04 

0.06 

0.08 

0.10 

0.12 

II    II    II    II 

0.00262 

0.03852 

0.03860 

-O.OOI8I+ 

0.00773 

0.07387 

0.07449 

-0.00468 

0.01351 

0.10605 

0.10806 

-0.00694 

0.01892 

0.13532 

0.13982 

-0.00788 

0.02317 

0.16182 

0.17029 

-0.00710 

0.02572 
0.18601 

0.19989 
-0.00452 

ry,   " 

-0.0002 

0 
-0.000003 

O.OOOli 

-0.0008 

0 
-0.000024 

0.0016 

-0.0018 

0 
-0.000081 

0.0036 

-0.0032 
0 

-0.000192 
0.0064 

-0.0050 

0 
-0.000375 

0.0100 

-0.0072 

0 
-0.000648 

0.0144 

"^■^aa  = 

-0.9336 

-0.8684 

-0.8073 

-0.7514 

-0.7014 

-0.6576 

'^"^ar'  = 

0.02369 

0.03209 

0.03480 

0.03419 

0.03150 

0.02752 

^~^a7"  = 

-0.U529 

-0.4094 

-0.3707 

-0.3371 

-0.3089 

-0.2862 

ar^y-y   = 

O.O292U 

0.05414 

0.07468 

0.09088 

0.10273 

0.11024 

0-2  z^  = 

0.02941+ 

0.05494 

0.07648 

0.09408 

0.10773 

0.11744 

"^"^W  = 

0.00063 

0.00251 

0.00566 

0.01005 

0.01571 

0.02262 

co-2z^^„  = 

lo. 02944 

0.05494 

0.07648 

0.09408 

0.10773 

0.11744 

^■^ifoa  = 

0.9631 

0.9233 

0.8837 

0.8455 

0.8091 

0.7750 

a)-2Zf^^,  = 

-0.02306 

-0.02958 

-0.02914 

-0.02414 

-0.01579 

-0.00490 

o^-^lf^y-  = 

0.48235 

0.4643 

0.4471 

0.4312 

0.4166 

0.4036 

2 

(a~        Ifyy      = 

0.0002 

0.0008 

0.0018 

0.0032 

0.0050 

0.0072 
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Figure  1.-    Elastic  suspension  of  a  Joukowsky  profile  with  plane  "spine." 
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Figure  2.-    Vectorial  lift  and  moment  coefficients  for  equal  beat  amplitudes. 
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Figure  3.-    Vectorial  lift  and  moment  coefficients  for  equal  beat  velocity. 
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Figure  4.-    Power  and  drag  coefficients  for  equal  beat  velocity. 

(1)   co-2zfaa;    (2)  a."2zfc.7';    (3)   a)"^fa7";    (4)    ao'Szf^r; 
(5)   (Jo"2wafl,;    (6)  Go"2wa7';    (V)   co'^Wa?";    (8)    03-2^77. 
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Figure  5.-    Degrees  of  efficiency  as  a  function  of  w    for  different  amplitude 

ratios    c. 
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Figures.-    Degrees  of  efficiency  as  a  function  of  the  phase   cp, 
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Figure  7.-    Test  arrangement. 
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Figure  8.-    Oscillation  diagrams. 
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